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We propose a “feature-scattering” mechanism to explain the cosmic microwave background cold 
spot seen from WMAP and Planck maps. If there are hidden features in the potential of multi-field 
inflation, the inflationary trajectory can be scattered by such features. The scattering is controlled 
by the amount of isocurvature fluctuations, and thus can be considered as a mechanism to convert 
isocurvature fluctuations into curvature fluctuations. This mechanism predicts localized cold spots 
(instead of hot ones) on the CMB. In addition, it may also bridge a connection between the cold 
spot and a dip on the CMB power spectrum at i ~ 20. 


I. INTRODUCTION 

In the recent years, observations of the cosmic mi¬ 
crowave background (CMB) radiation fluctuations by the 
Wilkinson Microwave Anisotropy Probe (WMAP) and 
Planck satellite have led to a precise measurement of 
temperature fluctuations on the sky from the largest 
scales down to arcmin scales mca- The temperature 
anisotropy is found to be highly Gaussian and “statisti¬ 
cally isotropic” in the sense that nearly all statistical pro¬ 
prieties of the temperature anisotropy can be described 
by the angular power spectrum Cj T [3]. However, it 
was found since WMAP 1-year data that there is a deep 
cold spot (AT ~ —120 K) in the southern Galactic hemi¬ 
sphere along the direction (l = 209°, b = —57°) with 
angular radius 0 ~ 10° mm, which is further confirmed 



FIG. 1: An illustration of the feature scattering mechanism. 
The classical inflationary trajectory is featureless. However, 
because of the quantum fluctuation of the isocurvature di¬ 
rection, hidden features may be encountered. At such an 
encounter, the inflationary trajectory gets scattered, the in- 
flaton losses kinetic energy and thus isocurvature fluctuation 
converts to curvature fluctuation. 


by Planck nominal mission data [3]. The cold spot is 
highly non-Gaussian in the sense that the probability of 
the cold spot existing in the statistical isotropic Gaussian 
universe is less than 0.1 percent [6]. 

Since then, the cold spot feature in the CMB map 
has invoked many observational and theoretical investi¬ 
gations. Initially, it was suggested that the unsubtracted 
foreground contamination might be responsible for the 
apparent non-Gaussian features m m, but later stud¬ 
ies mm showed that the significance of cold spot is not 
affected by Galactic residues in the region of the spot. 
It was also proposed that a spherically symmetric void 
with radius ~ 300 Mpc at redshift z = 1 can produce a 
large and deep CMB cold spot through the late-time in¬ 
tegrated Sachs-Wolfe effect (ISW) |10] (also known as the 
Rees-Sciama effect ED- Later studies PUCES] with the 
galaxy survey data El did find such a supervoid of size 
r ~ 195 Mpc with density contrast — —0.1 at redshift 
z = 0.16 align with the cold spot direction. However, 
more detailed following-up studies mm showed that 
the Rees-Sciama effect produced by such a void is sev¬ 
eral orders of magnitude lower than the linear ISW effect 
therefore is not able to account for the observed feature. 

The interesting non-Gaussian feature of cold-spot also 
invokes theorists to investigate the plausible explanation 
from the early universe. By considering various cosmo¬ 
logical defects in the early universe, Refs. El CEB] pro¬ 
posed that a cosmic “texture” (i.e. a concentration of 
stress-energy and a time-varying gravitational potential 
due to the symmetry-breaking phase transition) can gen¬ 
erate hot and cold spots on the last-scattering surface, 
with the fundamental symmetry-breaking scale found to 
be 0o ^ 10 15 GeV. However, by applying the Bayesian 
method to WMAP full-sky data, Ref. [19, did not find 
strong evidence of the texture model, neither completely 
rule out the possibility (at 95% confidence level). It was 
also proposed that cosmic bubble collision, predicted by 
eternal inflation theories, can induce the density pertur- 
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bation between our bubble and others, which can give 
arise to the localized features in the CMB [201 HU- But 
more detail data analysis [22] showed that the expected 
number of bubble collision is too few to account for the 
features in the CMB. Alternatively, a cold spot may fol¬ 
low from a different trajectory during multi-stream in¬ 
flation [23, [24] or modulated reheating after multi-field 
inflation [25] . 

The above physical or astronomical interpretations of 
cold spot either fail at some level, or require fine tuning 
or exotic scenarios of the early universe. Economically, 
some cosmologists would prefer to interpret the cold spot 
merely as a “3<r” statistical fluke. In this paper, we will 
provide a natural and physically plausible explanation of 
the cold spot, through multiple-field inflation. If the in¬ 
flationary trajectory is scattered by a feature hidden in 
the isocurvature direction, the inflaton loses some energy 
and thus inflation tends to be longer. Therefore, it is 
possible that only a small portion of the sky hits the fea¬ 
ture due to stochastic fluctuations, then that local patch 
of the universe experiences longer period of inflation and 
thus produces a cold spot. The mechanism is illustrated 
in Fig. [l] 

This paper is organized as follows. In Section [IT] 
we provide an explicit example of feature scattering, 
from massless isocurvature directions. Our predictions 
are compared with the measurement of the cold-spot in 
Planck's SMICA map. In Section [TTT[ we consider isocur¬ 
vature directions with mass m ~ H. We conclude in Sec¬ 
tion |IV] and discuss possible future directions. Through¬ 
out the paper, the unit M p i = 1/VSttG = 1 is used unless 
otherwise stated. 


II. COLD SPOT FROM FEATURE 
SCATTERING 


Such a potential can be derived from brane dynamics 
of string theory [26]. We will nevertheless not import 
anything other than this potential from string theory, 
but consider it as a phenomenological example. 

To fit the observations, we will take /i = 0.01 and Vo — 
4.89 x 10 -14 [27]. At the horizon crossing scale of the cold 
spot, the inflaton field is 0* ~ 0.153 for TV* — 4 (counting 
from the start of the observable stage of inflation), and 
the Hubble parameter is H* ~ 1.28 x 10 -7 . The feature 
SV is put at (j )o = 0*. 

In the isocurvature direction, the initial value of the 
field x* has no classical preferred value. We set (x*) = 0 
(one can always shift the xo parameter to satisfy this re¬ 
quirement if needed). The quantum fluctuation of x* (x) 
will be crucial in the analysis. Following the cosmic per¬ 
turbation theory, x* (x) has a scale invariant power spec¬ 
trum in Fourier space and we normalize the power spec¬ 
trum to be P x = [H/(2tt)\ 2 . In position space, X*( x ) 
behaves as a Gaussian random field. As we will show, 
the nonlinear mapping from x*( x ) to AT(x) generates 
the cold spot on the CMB. 

To explore the parameter space of SV, we will choose 
three sets of parameters as three illustrative examples of 
feature scattering. 

• Benchmark 1: negative short feature. We set 
A = -1.1 x 1O _5 VoM 4 /0*, a (f) = 10H*/(27r), g x = 
0.7H*/(27r) and xo = 4H*/(27r). The procedure of 
the calculation is also introduced in this subsection. 

• Benchmark 2: negative long feature. We set 

A = -1.4 x 10~ 5 VotJ, 4 /<J>i, = 50//*/(2tt), = 

0.7//*/(27r) and xo = 5.2//*/(2 tt). The differences 
from Benchmark 1 are discussed. 


A. Potential choice and method of calculation 


In this section, we provide an explicit example of fea¬ 
ture scattering to illustrate the physical mechanism. The 
generalization to other inflationary potentials should be 
straightforward. Consider inflation with two fields, both 
have standard kinetic term, and the potential is given by 


V 

SV 


v s M) + sv , 


A exp 


(4> - </>o) 2 

H 


(x - xo) 2 


( 1 ) 


The classical initial trajectory is chosen to be on di¬ 
rection, i.e. x = 0- Here V sv ((j)) is the slow-roll part of 
the potential. Our analysis is not sensitive in the form of 
V SY ((/)), as long as it successfully drives slow-roll inflation. 
For illustration purpose, we choose a small field potential 
for F sr (# 

VM=Vo(l~^j • ( 2 ) 


• Benchmark 3: positive feature. We set A = 1.0 x 
10~ 4 V 0 = 50Z/*/(2tt), a x = 0.5//*/(2tt) 
and xo = 3.5//*/(2-7r). The differences from Bench- 
marks 1 and 2 are discussed. 


Note that in the above benchmarks, the randomness 
in the potential needs not to be large. In fact, SV/V ~ 
10“ 5 is already significant enough to bring observable 
difference. This is because, the isocurvature fluctuation 
control whether the inflation trajectory hit the feature. 
And once the trajectory hits the feature, the change of 
curvature perturbation will be order ( ~ SV/V ~ 10 -5 . 
As a result, our mechanism is extremely sensitive to small 
features in the potential. 

The SN formalism 321134] will be used to investigate 
the cosmological perturbations from the feature scatter¬ 
ing. The SN formalism uses the observation that differ¬ 
ent Hubble patches during inflation can be approximated 
as different local FRW universes. Thus the cosmological 
observables can be calculated by exploring those local 
FRW universes. Especially, the curvature perturbation 
in the uniform energy density slice can be calculated by 
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FIG. 2: Left panel: 9n4> as a function of N. The red dashed line is for without feature scattering, and the blue solid line is with 
feature scattering. Note that the 4> field loses its kinetic energy very sharply at the 4-th e-fold (where the feature scattering is 
happening), and the kinetic energy slowly recovers after of order 1 e-fold. Right panel: The kinetic energy K(<f>)/H# = (<9jv</>) 2 
and K(x)/H* = ( &nX ) 2 - The color on the plot denotes e-folding number. One finds that about 0.1% of the 4> field kinetic 
energy transfers to x- Here we have taken x* = right before hitting the feature. 



FIG. 3: The temperature fluctuation as a function of x* — Xo- The left, middle and right panels are for parameters of 
Benchmark 1, 2, 3 respectively. 


C = SN, (3) 


where SN is the difference of e-folding number between 
an initially-flat slice and a finally-uniform energy density 
slice. The curvature perturbation ( is then converted to 
the CMB temperature anisotropy following the standard 
theory of CMB. In this paper we will use the Sachs-Wolfe 
approximation ST/T ~ — C/5 = —SN/5. Intuitively, the 


1 See the appendix of m for clarification of sign conventions. 


relation can be understood as: 

Inflaton scattered by feature 
-> smaller kinetic energy in </> direction 
-T larger SN 
-T later reheating 
-T less dilution of energy 
higher energy density 
-T deeper gravitational potential 
—> lower CMB temperature . (4) 

More precise relation between ST/T and ( can be studied 
by solving the Boltzmann equations at recombination. 

Note that SN has two sources: the quantum fluctua¬ 
tion of (p and the quantum fluctuation of x- Those inde¬ 
pendent sources can be studied independently, and they 
needs to be added together to make the total curvature 
fluctuation. In the numerical calculation, we will focus 
on investigating SN as a function of x but we will add 
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FIG. 4: Left column: The cold spot profile (we have rotated the map such that 0 — 0 corresponds to the center of the cold 
spot). The gray bars are 1000 simulations of fluctuations of the feature model. The green bars are simulation of isotropic 
Gaussian fluctuations (5000 simulated maps). Right column: five best-fitting cold spot profiles from the Benchmark models. 
The red dot is the actual cold spot from the Planck SMICA map for both panels. The three rows from top to bottom are the 
Benchmark models 1, 2, 3 respectively. Here the Sachs-Wolfe approximation is used when numerically fitting the bubble profile. 
The Sachs-Wolfe approximation works well only on large scales. Thus the fitting of the first two points with simulation may 
not be accurate enough. But one can see that the error bars of those two points from simulation are also large so that the 
inaccuracy should not affect the physical result significantly. 


back a contribution of Gaussian random field to account 
for the contribution from at the last map-making stage. 

The angular size of the cold spot is around 10 degrees. 
Converting to horizon crossing time of comoving wave 
number, this corresponds to i ~ 20. Note that in the 
observed CMB temperature power spectrum, there is in¬ 


deed a dip at i ~ 20. It will be interesting to do a com¬ 
bined analysis on both signals, from feature scattering. 

In the following subsections, we will explore different 
parameter space of Eq. ( 0 ) with three benchmarks, and 
calculate the e-folding number as functions of the posi¬ 
tion of isocurvature field SN = SN(x)- 
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B. Benchmark 1: Negative short feature 

In this subsection, we consider a negative short feature 
in the potential. We set A = —1.1 x 10 _5 VbM 4 /0*> = 

10tf*/(27r), <j x = 0.7tf*/(27r) and Xo = 4ff*/(27r). The 
size of the feature in the <p direction is relatively small 
(compared to the case of Benchmark 2). As a result, the 
a field does not have enough time to oscillate inside the 
negative feature. 

Before getting to the implications to the late universe, 
it is intuitive to check what was happening during infla¬ 
tion. 

As an example, we consider the X field value when 
hitting the feature to be X* — Xo = H*/(2it). The time 
evolution of is plotted on the left panel of Fig. |2j 


One can find from the plot that, at TV* ~ 4 (where feature 
scattering happens) 


dc/) 

d N 


4 x 1(T 4 . 


( 5 ) 


When the inflationary trajectory hits the feature, there 
is a change in the velocity of the inflaton </>. At first, 
\d(/)/dN\ increase because </> falls into a potential well. 
However, two other effects immediately follow: First, the 
obtained energy is returned because (j) has to climb out of 
the potential well; second, the potential well can be con¬ 
sidered as a scattering center, such that a similar amount 
of energy is transferred to the x direction (as evidenced 
in the right panel of Fig.j2|. As a result, there is a overall 
loss of velocity of the inflaton, with the value read from 
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the plot 


dc/) 

d N 


- -4 x 10- 7 


- 10“ 3 


dc/) 

d N 


( 6 ) 


In words, the inflaton field has lost 0.1% of its kinetic 
energy because of hitting the feature in the potential 2 . 
Note that the inflaton always loses its kinetic energy be¬ 
cause of scattering off a feature. Thus a cold spot instead 
of a hot spot is predicted 3 . 

Note that the loss of kinetic energy takes of order one 
e-fold to recover, because this is the time scale to reach 
the inflationary attractor solution. To be more precise, 
as shown in the left panel of Fig.[2j it takes about roughly 
0.2 e-fold for the inflaton to recover its kinetic energy. As 
a result 4 , 


SN = 



0.2 x 

d 0 
dN 


&4> 

dN 

+ s 

dcf) 

dN 


2 x 10“ 4 . (7) 


From the Sachs-Wolfe approximation, the temperature 
fluctuation 5T/T ~ — £/5 = —8N/5. Thus we get 5T ~ 
— 100/iK. 

Carrying on the above analysis for general values of 
X* — Xo, numerically the temperature fluctuation as a 
function of x* — Xo is plotted in the left panel of Fig. [3] 
We then simulate 1000 maps and rotate those maps such 
that they center at the cold spot. In the left panel of the 
top row of Fig. [4j we bin the pixels of the map as a func¬ 
tion of 0, and then plot the 0.5, 0.16 and 0.84 quantiles 
(corresponding to the central value and la uncertainty if 
the probability distribution were Gaussian) of the binned 
pixel temperature. In the right panel of the top row of 
Fig-U we plot the top five best-fitting cold spot profiles. 
Three best-fitting examples are given in the top row of 

Fig.0 


C. Benchmark 2: Negative long feature 

In this subsection we consider a negative long feature. 
Inside a “long” feature, there is enough time for the x 


2 Here we do not need to consider the change of potential energy 
before/after the feature scattering, because we only need to com¬ 
pare the change of energy before and after a very sharp scattering 
process, and the change of potential energy is negligible. 

3 When the potential is positive, the inflaton actually gains kinetic 
energy at the very beginning, when it falls into the potential well. 
However, the gained energy is returned after a time interval that 
is much smaller than Hubble time. Thus the effect of gaining 
the kinetic energy is negligible when we integrate to obtain the 
e-folding number. On the other hand, the loss of kinetic energy 
through feature scattering needs order of 1 e-fold to recover, and 
thus it is the dominate effect. 

4 The 5N here is from the quantum fluctuation of %. The part 
from (j> will be added at the map-making stage. 


field to oscillate. We set A = —1.4 x 10 -5 Vb///</> 4 , a$ = 
50H*/(27r), a x = 0.7H*/(27r) and xo = 5.2H*/(27r). 
Note that the feature size in the <fi direction is 5 times 
longer than that of Benchmark 1. 

Finer structures develop in the case of long feature. 
If we have chosen an initial value X* — Xo = H*/( 27r), 
we find similar behavior to Fig. [2j However, for some 
other initial values, for example, X* — Xo = 1.4i7*/(27r), 
the x field oscillates inside the feature. The situation 
is illustrated in Fig. [6j As a result, some amount of x 
kinetic energy is returned to </>. For some special values, 
almost all kinetic energy is returned. Thus oscillatory 
patterns develop, as shown in the middle panel of Fig. [3] 

The oscillations in ST(x) leads to ring objects in 5T(x). 
The cold spot in this parameter space is not completely 
cold, but instead has nested rings of cold and hot pat¬ 
terns. There is no evidence of such patterns in the actual 
CMB cold spot. There has been a debate if there are 
other ring patterns in the sky |36j . The Benchmark 2 pa¬ 
rameters provide an explanation for such rings (though 
the shape is not perfectly spherical). However, it is shown 
that the appearance of such rings is due to an inappro¬ 
priately chosen power spectrum, and thus not actually in 
the sky [37 . Thus we will not tune the parameters to fit 
such features here. 

The general and best-fitting cold-spot profiles, and 
sample spots are plotted in the middle row of Figs. [4] 
and[5j respectively. 


D. Benchmark 3: Positive feature 


In this subsection we consider positive feature scatter¬ 
ing. We set A = 1.0 x 10 _4 VbM 4 /^*> = 50H*/(27r), 

a x — 0.5i7*/(27r) and xo = 3.5H*/(27r). With positive 
potential, the x direction has a run-away solution and 
can never be bounded. So there is no oscillation in this 
case. 

The ST(x) dependence is not as sharp as the previous 
two benchmarks. Thus the cold spot does not have as 
clear a boundary as those previous cases. Because of 
the not-so-sharp ST(x) dependence, to get a cold enough 
center of the spot, we have to increase A. As a result, the 
cold spots generated by Benchmark 3 typically (though 
not always) have very cold tinny cores. 

The general and best-fitting cold-spot profiles, and 
sample spots are plotted in the bottom row of Figs. [4] 
and[5j respectively. 

Considering a very sharp and cold core is formed for 
positive features, it is unlikely that such features can fit 
the realistic cold spot. But considering that the sharp¬ 
ness of the cold core depends on model parameters. The 
cold core may appear at £ > 2000 scales and free stream¬ 
ing may erase the sharpness. 
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FIG. 6: The same as Fig. [2] but for Benchmark model 2. Note that the feature in the inflationary potential is long enough to 
hold a few oscillations in this case. 



6 [degree] 



6 [degree] 


FIG. 7: The same as Fig. [4] but for massive isocurvature model m x = 0.9 H*. 


III. MASSIVE ISOCURVATURE DIRECTIONS 

The feature scattering not only applies to massless 
isocurvature directions, but also to marginally massive 
ones. As long as the mass of the isocurvature direction is 
not much greater than i7, the energy scale for inflation¬ 
ary fluctuations, the isocurvature direction can still be 
excited to explore its field space. During inflation, fields 
with m x ~ H arise naturally [281431] . Thus it is worth¬ 
while to investigate feature scattering of such marginally 
massive fields, with an additional mass term in the po¬ 
tential 

Vnass(x) = \ m2 X X 2 ■ (8) 

Note that the process of feature scattering happens on 
a time scale much shorter than Hubble. Thus as long as 
m x is not much great than H*, m x does not enter the 
dynamics of feature scattering. However, the spectral 
index of a massive field is significantly different from a 


scale invariant spectrum, with spectral index 

n *- 1 = Re (^i-S) ■ (9) 

Thus the mass of the field controls how many e-folds can 
the isocurvature direction randomly travel. For example, 
when m = 0.9H*, the x field has a spectrum n x — 1 = 0.3. 

In the numerical example, we set m = 0.9i7*, A = 
-1.1 X 10 ~ 5 V 0 n 4 /(t)i, a 4 , = 10iJ*/(27r), a x = 0.7ff*/(27r) 
and xo = 4.4ff*/(27r). In other words, we have only 
tuned xo to be slightly greater, and consider a massive x 
field. Otherwise the parameters coincide with the Bench¬ 
mark 1 of massless case. 

The general and best-fitting cold-spot profiles, and 
sample spots are plotted in Figs. [7] and [8j respectively. 

As one can see from Fig. [8] there are more small-scale 
structures compared to the case of massless Benchmark 

1 . 
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IV. CONCLUSION AND DISCUSSION 

To conclude, we propose a scattering mechanism of in¬ 
flaton due to the features in the inflation potential, dur¬ 
ing which the isocurvature fluctuations are converted into 
curvature fluctuations. The curvature fluctuations direc¬ 
tion loses kinetic energy due to the scatters in isocurva¬ 
ture direction, therefore the number of e-folds becomes 
larger in some region of the universe. We find that the 
cold spot can be well explained by such a mechanism, 
and the spot profile reasonably fits the CMB cold spot 
without fine tuning of the inflationary parameters. Be¬ 
fore ending up the paper, we would like to mention a few 
future directions: 

• Beyond the Sachs-Wolfe approximation. Currently 
we did not use the full Boltzmann code to calcu¬ 
late the CMB transfer function. The Sachs-Wolfe 
approximation works well for exploring the coarse¬ 
grained cold-spot profile but is not enough for prob¬ 
ing the fine structures inside the spot. One can use 
the full radiation transfer function to carry on a 
more detailed analysis in the future. 

• Verifying the approximations of the SN formal¬ 
ism. The SN formalism assumes that the horizon¬ 
crossing amplitude of the 0 and x quantum fluctu¬ 
ations are Gaussian and uncorrelated. In the case 
of feature scattering, this is not rigorously true be¬ 
cause of the turning of trajectory. However, note 
that the transfer of the inflaton kinetic energy is 
tinny (0.1% in the studied example), the correction 
from sub-horizon physics should be suppressed by 
this small fraction. 

Having that said, it remains interesting to see if the 
calculation from SN formalism can be verified by 
first principle calculation. However, such a calcula¬ 
tion is challenging. The first principle calculation 
of cosmological perturbations is known as the in-in 
formalism. But in the in-in formalism, the primary 
calculatables are the correlation functions. In other 


words, one starts from the Gaussian two point cor¬ 
relation function and studies small departure from 
that. But the cold spot is highly non-Gaussian and 
localized object, and thus is not easily captured by 
the correlation functions from the in-in formalism. 

• Extra species/symmetry point (ESP) [38] 39] as 
scattering centers. In our current examples, the 
inflaton kinetic energy is lost into the collective mo¬ 
tion of the isocurvature direction. It may also be 
possible that the kinetic energy of the inflaton is 
lost into hidden ESPs in the isocurvature direction. 
It is interesting to explore such possibilities. 

• Connection between the feature scattering regime 
and the multi-stream regime: By carefully arrang¬ 
ing the position of the feature, it is possible that 
different part of the universe follows different clas¬ 
sical trajectories, separated by temporary domain 
walls. This is known as the multi-stream inflation 
[23, .24;, which is also a possible explanation of the 
cold spot. Feature scattering and multi-stream cor¬ 
responds to different regimes of the multi-field pa¬ 
rameter space. If the features on the inflationary 
potential are random, we expect the feature scatter¬ 
ing to be more typical (whereas highly random fea¬ 
tures in multi-stream inflation cause disasters and 
provide a constraint on multi-field inflation [40]). 

• Connection between the cold spot and the dip of 
CMB temperature power spectrum at £ ~ 20. We 
argue they may come from the same origin - fea¬ 
ture scattering during inflation. With our current 
toy potential given at Eq. (Q), we are unable to 
explain this power deficit at the same time. How¬ 
ever, a more general shaped potential may be able 
to explain both the deficit and the cold spot. 

• String theory model building. It is believed that 
string theory has a landscape of complicated vac¬ 
uum structures [41] [42]. It would be interesting to 
build string landscape models for feature scatter¬ 
ing. Especially, our mechanism assumes that 
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and a x have value of order H. This assumption is 
made at the phenomenological level in this paper. 
Nevertheless, it is interesting to search for theoret¬ 
ical motivations for such a potential. For inflation 
to work, the inflationary perturbations, in the sense 
of effective field theory, should have enough range 
of validity. Especially, the UV cutoff of the theory 
should be at least of order H. In the marginal case 
where the UV cutoff of the effective field theory is 
just of order H , new features are expected when 
the field rolls a distance of order H . In addition, 
inspired by quasi-single field inflation |28j [29] , the 
mass parameters during inflation should be natu¬ 
rally of order H. Thus once there are many mass 
parameters in the landscape, such as some mass 
matrix, variations will be introduced in the poten¬ 
tial with field range of order H. We leave detailed 
investigation of those possibilities to future works. 



FIG. 9: Relation between xo and cr x , in order to keep the 
cold spot(s) covering 0.76% of the sky. The xo axis is plotted 
with unit c — H*/ ( 27 t). 


• Careful model comparison between our feature 
scattering mechanism and the ACDM base model. 
We have 5 parameters in our illustrative model. 
The parameters are chosen to be the convenient 
ones for inflation model building, but not opti¬ 
mized to fit the data as a phenomenological model. 
It remains interesting to propose a more econom¬ 
ical model with fewer number of parameters, and 
perform model comparison between our mechanism 
and standard inflation model. 


To be more explicit, the relation between a x and xo is 

r i (x-xo ) 2 

/ — l= e 2c2 dx ^ 0.0076 , (Al) 

J-* x cv27t 

where c = y/ (x*)- Here we have used because there 
are uncertainties from both different realizations of ran¬ 
dom variables, and the detailed profile of the inflationary 
potential (considering the feature does not sh arply start 
at — a x and sharply end at a x ). Equation ( |A1| ) does 
not have an analytic solution. The numerical solution is 
plotted in Fig. [9] and can be fitted to high precision by 
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' Xo ~ 2.43c, 


< -Xo - 2.43c, 

. [exp (^^) 


0.99 


if xo > 2.8c 

if xo < -2.8c 

if — 2.8c < xo < 2.8c 

(A2) 


This relation is plotted and referred to as “0.76% fit func¬ 
tion” in Fig. [9j In the num erical examples, a x and xo 
are chosen around Eq. (A2), and slightly adjusted con¬ 
sidering the detailed profile of the potential. 


Appendix A: Generation of Gaussian random fields 

The fluctuation of the isocurvature direction follows 
Gaussian statistics. In this appendix, we provide some 
details of statistics. 

The cold spot has an area of about 0.76% of the full 
CMB sky. To have such a fraction of sky covered, there 
are two possibilities: Either the feature is very narrow 
(a x < i7*), or xo is about 2.5cr away from a x . Practi¬ 
cally, the former tends to produce a number of small cold 
spots and the latter tends to produce one or a few large 
cold spot (which is the case of our current interest). 


Appendix B: Selection of spots 

In the numerical simulation of the CMB maps, the cold 
spot is selected with the following pipeline: 

• The simulated sky is a 200x200 square degrees. 
This approximately resembles the whole CMB 
sphere. Note that we will eventually focus on the 
spot, which is around 10x10 square degrees, thus 
the flat sky approximation can be applied here for 
technical simplicity. 

• For each map, we generate a mask, where the un¬ 
masked pixels have temperature ST < — 100/iK. 
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FIG. 10: Samples of rare fluctuations out of a Gaussian random field. The upper and lower panels show samples of 2.5 a and 
3 a rare fluctuations respectively. 



• The mask is smoothed with a Gaussian filter with 
a width of 2 degrees. After Gaussian smoothing, 
the very nearby ones are connected. 

• The largest connected (after smoothing) spot on 
the mask is selected. 


• We crop the real image according to the spot po¬ 
sition of the mask, with a radius of 30 degree. We 
drop the images if the spot is too close to the 
boundary of the whole image. 

Some examples of spot shapes are illustrated in Fig.[l0| 
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